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. We calculate the entropy of adiabatic perturbations associated with a truncation of the 

O^l ' hierarchy of Green functions at the first non trivial level, i.e. in a self-consistent Gaussian 

approximation. We give the equation governing the entropy growth and discuss its phe- 
' nomenology. It is parameterized by two model-dependent kernels. We then examine two 

particular inflationary models, one with isocurvature perturbations, the other with correc- 
tions due to loops of matter fields. In the first model the entropy grows rapidely, while in 
the second the state remains pure (at one loop). 
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I. INTRODUCTION 
Oh' 

, This is the second part of a series of notes on the decoherence and entropy of primordial 

fluctuations predicted by inflation. The first part, subsequently called I was devoted to 
the operational formulation of the notion of decoherence of metric perturbations. Here we 
^ , turn to the calculation of the entropy when truncating the hierarchy of Green functions 

I at the first non-trivial level, first in general and then in two particular models. As an 

■ introduction, we briefly review the present state-of-the-art concerning the entropy growth 

during inflation. 
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I Up to now, the time dependence of the entropy as well as the main source(s) of entropy 

■ remain undetermined. Part of the difficulties arises from the definition of the reduced 

. density matrix. Concerning the methodology, previous studies either work with a master 



X 



equation for the reduced density matrix [3, y] , or directly calculate the reduced density 
matrix [1, [E], or use analogies with other quantum mechanical situations 0|. 

The main problem with the first two approaches comes from the actual calculation of 
5^ I the trace. While this poses no difficulty for Gaussian models as in [1, [H], it is hamperred 

by the infinities ubiquitous in interacting Quantum Field Theories. In previous attempts 
where this calculation was undertaken [Sj], unjustified assumptions were made and lead 
to contradictory results. We argued in paper I that the only way to deal properly with these 
infinities is to work with Green functions. The approach of [g] is different as it assumes a 
general Gaussian and Markovian Ansatz for the master equation. The limitation of this 
approach is that it does not make predictions. Indeed, first the parameters of the master 
equation are undetermined unless they are calculated in field theoretical settings. Second, 
the Markovian hypothesis is unjustified for superhorizon perturbations at any epoch, as 
well as for perturbations of any scale during inflation. 

In this paper, we present an approach which combines the advantages of the previous 
ones, but without their shortcommings. We work with Green functions from the onset. 
Infinities are handled by the standard renormalization techniques, and the general evolu- 
tion equation for the entropy follows straightforwardly from the knowledge of these Green 
functions. This emphasis on Green functions is not new [3], but to our knowledge it has 
not been applied to the calculation of the entropy of primordial fiuctuations in realistic 
models of infiation. 
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After presenting the settings in Sec. [Til the evolution equation of the entropy of the 
reduced density matrix is derived in Sec. IIIIl followed by a general discussion of this 
equation. We then examine two specific models. A model of two field inflation is presented 
in Sec. IIVI It confirms and generalizes the analysis of [H]. In particular we identify the 
conditions under which the entropy grows with the number of efolds. In Sec. |V] we 
compute the entropy associated to one-loop corrections from matter fields in a class of 
semi-realistic theories analyzed by Weinberg in [sl. In these models we show that no 
significant decoherence occurs during inflation. 



II. THE GAUSSIAN APPROXIMATION 

We work in the coordinate system in which the inflaton fleld is homogeneous on each 
space-like hypersurface, 

v{t,^) = Mt) ■ (1) 

We focus on the curvature perturbations ^. In the linear approximation, their evolution 
is described by the action 



dtd^xa^lC^ - ^{VCY] , (2) 



5 = ^/.«aa(c^-J,'-= 

where a is the scale factor and e = —H/H^ is the first slow roll parameter, see for a 
comprehensive derivation. A dot stands for the derivation w.r.t. the cosmological time t. 
The Fourier mode labeled with the conserved comoving wave vector q obeys 

C, + ^^C, + fjC, = o, (3) 

where q is the norm of the vector q. 

The free vacuum is taken to be the Bunch-Davis vacuum, defined from the positive 
frequency solutions of ([3]) for infinite physical momentum, 

(i9.-g)(aV^Cr) -0, 

r is the conformal coordinate time dt = adr. Using this solution, the mode operator can 
be decomposed as 

Cq(t)=a^°CrW+«^q^C*W^ (5) 

where the destruction operator annihilates the Bunch-Davis vacuum. 

We shall work in a quasi-de Sitter approximation. It can be shown that the results are 
still valid in slow roll inflation. In this approximation, the in mode is 

= C°(l + igr)e-^^. (6) 
The constant is fixed by the equal time commutator [^q, vr^,] = i5^{(\ — q'), where 



is the Fourier component of the momentum conjugate to Q. One finds 

47rG H"^ 
2^ 
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In the hnearized description based on Eq. ^ , the state of ^ stays Gaussian and its prop- 
erties are characterized by a single function, the power spectrum oc In interacting 
field theories, the state is characterized by the full hierarchy of connected Green func- 
tions. Such knowledge is out-of-reach, so that in practice one resorts to a (self-consistent) 
truncation of this hierarchy. This coarse graining defines a reduced density matrix with a 
non vanishing entropy, see Sec. Ill in paper I. The first non trivial level of truncation is 
defined by setting to zero all connected Green functions of order > 3. The reduced density 
matrice pred is therefore Gaussian. Since the state is still homogeneous, each two mode 
sector (Cq,C-q) can be analyzed separately. Therefore Pq^'lq is fully characterized [H by 
the c-number function G{T,T';q) given by 

G(T,r';g)53(q_q') ^ ilY {Cq(r), C-q'(r')}) • (9) 

Since the different sectors do not mix, we shall work at fixed q and no longer write the 
trivial (5^(q — q') coming from the plane wave normalization. 

To calculate the entropy carried by Pq^q it is convenient to recast the information 
contained in G(r, r'; q) into the covariance matrix C defined by 

As shown in Appendix I A 11 in the Gaussian approximation it is always possible to make 
a canonical transformation (Cj^r') i— > (C)^) where the new momentum is related to C as 
in d?]). This canonical transformation leaves invariant the entropy [l[. Using Eq. ([7]), the 
three moments V are then related to G by 

V(;{q,t)=G{t,t;q), (11a) 
rcn{q,t) = ^dt>G{t,t';q)\t=t> , (lib) 

/ 3 \ 2 

K(g,t)=f|^j dtdt>Git,t';q)\t=t'. (11c) 

For the linear perturbations in the Bunch-Davies vacuum, since G{T,T';q) = 
^{Cq^i'^) Cc^*{'^')) where is given in Eq. ([6]), the (unperturbed) covariances are 

" =|C°l'(l + x2) , (12a) 





= K{r) 






47rG 









\At:G) 

= ^vic°i^(|^)', (12c) 



where we used the relation r ~ —1/aH and introduced 



..-i^.e-. (13) 

= — Inx is the number of e-folds with respect to horizon exit. 

The entropy of Pq^-q is related to the determinant of the covariance matrix by 

S = 2 [(n + 1) ln(n + 1) - n ln(n)] , (14) 

where the parameter n is defined by 

1' 



n+- ^dei{C)=V^V.-Vl^. (15) 
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The prefactor 2 in (jl4p accounts for the fact that p^(^_^^ is the state of two modes. The 
pure states, S = n = 0, correspond to det(C) = 1/4. At any time, one verifies that 
Eqs. (jl2p exactly give this constant value of the determinant. Notice also that, had we 
neglected the terms 0{q^ /a?!!'^) in (fT^ ). we would have instead obtained det(C) = 
which makes no sense in quantum mechanical settings, since the inequality det(C) > 1/4 
is nothing but the Heisenberg uncertainty relations. 

We will also employ another useful parameterization of the covariance matrix 

nc(t) = Tr(p-^„q(t)af t«BD^ ^ ^^g) 

cc(t) ^ Tr a^D^ , (17) 

riQ is real while is complex. We recall that these parameters depend on the choice of 
canonical variables, that is on the choice of (aq,aq), see [H for more details. So does the 
parameter 5^ defined in the next equation. This parameter characterizes the correlations 
between the modes q and — q in the state pred^ 

|c^|2 = n^(n^ + l-(^^), 0<5^<n^ + l. (18) 

In this parameterization, the determinant of C is 

det(C)= (n,- + iy-|ccP = ^ + nc5c. (19) 

The entropy is a monotonously growing function of 5^. This parameter quantifies the 
residual coherence in the state: the state is pure when (5^ = and thermal when 5^ = n^+1. 
The threshold value 5^ = 1 separates entangled states from decohered states which cannot 
be distinguished from statistical ensembles [ij. At the threshold value, in the radiation 
dominated era, the entropy is 

5sep = In ~ - In (x^^^) = 4 iVcnd , (20) 

where A'end = — In x^nd is the number of efolds from horizon exit to the end of inflation. 
We have used the linearized treatment to estimate n™'^ (i.e. = l/4Xg^j) [13]. The 
maximal (thermal) entropy (per two-modes) is 5max = 2S'sep — 21nncnd — ^^cnA- 



III. EVOLUTION OF THE ENTROPY 
A. Equation of evolution 

As seen from Eqs. ()14p and ()15p . is a monotonically growing function of det(C). In 
particular, for 1, the entropy is simply given by 5 = ln(n(5) = In(detC). Hence, in 

this regime one has 

n(5 » 1 , ^=J^ • (21) 

The evolution of det(C) follows from that of the three expectation values of Eqs. ()lip 
which are all determined by the anticommutator G. Therefore the time dependence of 
det(C) is governed by the equation obeyed by G. The latter is a linear second order 
integro-differential equation. It is derived in Appendix lA 1[ But, as explained in lA 21 it 
is simpler to exploit the Gaussianity and derive G from an equivalent quantum Langevin 
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equation. That is, G is the anticommutator of the operator Cq(^) which verifies the effective 
equation 



Cq + 



d {Hah)) .^ ^ g2 ^^^^^ ^ rt ^^^^^ ^^^^^^^^^ ^ ^^^^^ ^^^^ 
^ J —oo 



dt 



In the Gaussian approximation, the effects of the interactions are summarized in two 
kernels, namely the "dissipation" kernel Dq and the so-called noise kernel 

Nq{t,t') = ^{{Ut),C-^{t')}). (23) 

N and D appear in the effective action of the curvature perturbation as the real and 
imaginary parts of the renormalized self-energy dSM]). For the differential part of ([22]) 
we assumed the same structure as in the free mode equation ([3|). Equation (j22p could 
be generalized to the case where the renormalized frequency differs from g^/a^. This 
would only affect the retarded Green function and G but not the form of the equation 
(f28l) governing det(C). However our result does not generalize to time dependent "wave 
function renormalization" since, should they occur, they would invalidate the use of Eqs. 
dMD and (125]). 

We now proceed assuming that the kernels D and N are known, and that G and its 
time derivatives possess well-defined coincidence point limits. We derive an equation for 
the covariances of For this we need the identities 

f = '^«' 

^ = P« + i ({C,(*). C-,(())> ■ (25) 

The equations can be easily derived in the Heisenberg picture, since time derivation and 
taking the expectation value, being both linear operations, commute. These equations are 
easily verified for ([12]) and Inserting (p2]) into (p5]) we get 

+ j"jt'Dqit,t') \ ({Cq(t), C-q(t')}) - C-q(t)}) • (26b) 

The three terms in (j26b ) are local in time and do not contribute to the time evolution of 
det(C), as will be seen in Eq. (128]) . Instead, the terms in ([26b ) are non local and absent 
in the free evolution. Without them, the time evolution of V^^ is entirely determined by 



that of V^^ and Vc^Ci as shown in Eq. (I25p . 

To compute the growth in entropy, we also need the time derivative of 7-"^!^, obtained 
from using again ([22^ . 



dt ~ ^i^q' ^ dt a2^C 

Jt'Dg{t,t') ({Cq(t), C-q(t')}) + {{Ut), C-q(t)}) • (27) 
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Using the equations d?]), (p6|) and ([27|) . we find 

+ r (it' [Ng{t,t')dtG,,t{t,t') - Dg{t,t')dtG{t,t')] \ , (28) 

where G^et is the retarded propagator associated with (j22p . The prefactor (a^e/AirG)^ 
comes from the relation ([7]). We will verify on several examples in Sec. HVl that this 
prefactor yields a rapid growth in (aH/q)^. 

Equation ([28j) governs the rate of change of the entropy defined by reducing the state 
as explained in Section [III It is generic and exact. The explicit expressions of the noise 
and dissipation kernels are of course model dependent. The important point is that they 
can be calculated using standard QFT techniques. 

This equation is equivalent to a master equation, but has several advantages over 
it. The master equation is the evolution equation for the reduced density matrix. It is 
obtained by calculating the partial trace on the r.h.s. of the Heisenberg equation 

i|p-% = IVq,^q([//,p]), (29) 

where H are p are the Hamiltonian and density matrix of the entire system. The trace is 
performed on all field configurations except (^q. The difficulty with (I29p (or the equivalent 
equation for its Wigner representation) is to actually calculate the trace. This is relatively 
straightforward in non relativistic quantum mechanics where the master equation has 



indeed proven to be a fruitfull approach, see [Ul] and [Ij, UM recent reviews. For 
quantum fields however, the first nontrivial contributions to the r.h.s. of (j29p arise from 
loop corrections. Therefore, the calculation can only be reliably done at the level of Green 
functions. Previous attempts to calculate directly (p9]) involved dubious assumptions and 
approximations which have lead to contradictory results, compare e.g. aaa. Moreover, 
Eq. (|29p does not bring us anywhere nearer to the growth of entropy since one still has 
to calculate the covariances. In contrast, Eq. ()28p involves the renormalized expectations 
values of the relevant observables. 

In the remainder of this section we adopt a phenomenological approach and discuss 
how S depends on the properties of N and D. We start with two general remarks. 

As a consistency check we verify that det(C) is constant in an equilibrium state and 
in a static universe. In that case, N and D and G and the commutator G[ j are related by 
the fluctuation-dissipation (or KMS) relations (for each q), 

iV(a;) = coth(^)z)(^), (30a) 
G(a;) = coth(^)aet(^), (30b) 



where T is the temperature. Inserting these identities into the r.h.s. of ()28p yields, as 
expected, det(C) = cte through a detailed balance. 

We also note that out of equilibrium, the sign on the r.h.s. of ()28p is a priori undeter- 
mined. However, one expects on physical grounds that the entropy averaged over sufficient 
time scales does not decrease.^ 



This is indeed what happens, for instance, during the early stages of the far-from equihbrium evolution 
of a quantum field [ij- Then the value of det(C) oscillates around a monotonously growing mean value. 
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We now discuss two limiting cases where N and D are approximately local in time or 
not. We then return to the general discussion of (p8]) in light of these results. 



B. Markov approximation 

This is by definition the regime where the correlation time of ^ is short compared to the 
characteristic time scale(s) of evolution of the cosmological perturbations. In this limit, 
both N and D are local in time. Let us assume that the environment resembles a thermal 
bath (contrary to the discussion at the end of the previous section we do not assume here 
that C is at equilibrium). We write the Ansatz 

Nit, t') = Mimt - 1') , D{t, t') = m^^it - 1') , (31) 

where M{t) and I'(t) are slow functions of time compared to the correlation time scale 
of the environment, caused for instance by the adiabatic expansion of the universe. Ac- 
cording to Eq. ([22]) dissipation corresponds to P > 0. Using the identities Gret{t,t) = 0, 
dtGret\t'=t = 1, and the definitions (fTT]) . Eq. (f28|) yields 



ddetjC) 
Jt 



/ 3 \ 2 

-2V{t)det{C) + 2(^] M{t)V^. (32) 



It is interesting that in this limit, the dissipation kernel acts on the determinant of C, 
whereas the noise kernel acts only on the power spectrum times the factor a^. To continue 
the discussion, we simplify Eq. (|32p by dropping the term proportional to T). This is a 
bad approximation to describe the approach to equilibrium, since (p2]l would not lead to 
det(C) — > cte. Far from equilibrium, it leads to an upper bound on the entropy (since P 
and det((7) are both positive). Then, for n 3> 1, the rate of growth of the entropy per 
efold is 

dS_ ^ a^{t)e^H~Wit)Vcit) 

The relevance of the Markov limit for cosmological perturbations is on the contrary dis- 
putable. The typical context where the Markovian limit emerges is that of an environment 



in a thermal state at high temperature [Ul]. This is irrelevant for inflation (both single and 
multi-fields) since the perturbations are in a squeezed state. This is also unlikely to be a 
realistic model after horizon reentry but before decoupling, because the scales entering the 
horizon are not in thermal equilibrium (since we observe acoustic peaks). We will come 
back to this point in the conclusions. 



C. Non Markovian regime 

At low and vanishing temperatures the Markovian approximation fails and all the terms 
in the r.h.s. of (|28|) are a priori of equal importance since D ^ N . Most likely this is 
the case relevant for (single field) inflation since all energy densities associated to particle 
excitations are redshifted away (with the exception of multifield inflation treated in Sec. 

When N and D are of the same order, we see on (j28p the possibility of a partial cancel- 
lation between the terms proportional to Vc_ and "P^^. It is perhaps easier to understand 
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this from the decomposition of the covariance matrix into the sum of the free contributions 
Eqs. ()12p and a remainder M 

C = Ci,,, + M. (34) 

This separation follows naturally from the perturbative calculation of G in (jlip where 
M regroups the corrections due to interactions, see Appendix I A 1[ In the representation 
(jl6til7p . the decomposition ([Mj) corresponds to 

= "-free + « , C = Cfree + " X • (35) 

where jcfreeP = '^free {nfree + 1)? where a is real and x is a complex number with norm 
smaller than 1. The quantity a is essentially governed by the power of the fluctuations of 
the environment, while x accounts for a possible squeezing of the state of the environment 
(which implies in particular N{t,t') / N'{t)5{t — t')). The corresponding value of 5 as 
defined in (fTB|) is 

5ia,x)=2a(l + \^) \l - Re (xe~'-'^^'^-^A] + ^{1 - \x\') + O ( ^) .(36) 

This equation clearly shows the competition between the contribution a, that tends to 
increase the entropy, and the contribution of x when the latter has a phase similar to 
arg(cfrcc)- In particular, when arg(x) = arg(cfree), the first term in (f36l) vanishes, 6 is 
quadratic in a, and will thus still obey (5 <C 1. Equation (j36p shows that due to subtle 
interference effects, decoherence may be largely suppressed when the perturbed state is 
squeezed as the unperturbed one. We will have an example of this phenomenon in Sec. 
IVl where the noise is due to matter fields in the interacting vacuum. Given the fine tuned 
nature of the condition arg(x) = arg(cfree)5 it is likely that this is the only case where it 
occurs. We notice finally that a similar mechanism has long been envisaged to circumvent 
the fundamental limit of quantum noise (i.e. the shot noise and pressure noise of the laser) 
in interferometric detectors of gravitational waves ^. 

D. Constant rate of entropy growth 

One can easily show that the entropy rate depends crucially on how fast the solution 
of (I22p asymptotes to a constant (when it does) outside the horizon. To do this, we use 
the equivalence between the quantum and stochastic versions of Eq. (j22p as explained in 
Appendix I A 2\ and we consider the Ansatz 

Ut) = C° (l + ax + /3y + 7 Y + 0(x^)) . (37) 

written with the variable x = q/aH. The coefficients a,... of the expansion are the sum 
of two terms, a = a + a: a constant term a = ((a)), corresponding to the positive 



^ In the third generation of those detectors, the numerous sources of noise could possibly be reduced 
enough so that the dominant one would be the quantum noise of the laser monitoring the position of 
the mirrors, see for instance [l^. In conventional interferometers, the laser produces two types of noise, 
the shot noise proportional to the intensity and the radiation pressure noise inversely proportional to 
the intensity. The total quantum noise can therefore not be arbitrarily small. This minimal value is 
called the Standard Quantum Limit. Noticing that light is squeezed during its travel in the arms of the 
interferometer, Unruh [3l proposed that judiciously squeezing the light send in the so-called dark port 
of the interferometer could beat the Standard Quantum Limit. 
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frequency solution ([6]) of the free mode equation, given by a = 0, /? = 1/2 and 7 = i/3; 
and a stochastic component d, ... which depend hnearly on the stochastic source and the 
dissipation kernel D, so that ((d)) = 0, .... Of course this Ansatz does not cover all the 
cases possible. We will consider another behaviour below. 
From ([37|) we deduce 

det(C) = ^ (^^^ [co + 2cix + C2X^ + c^x^ + 0(x4)] (38) 
where the coefficients are 

Co = ((!«!')) (39a) 

ci = ((Re(a*/3))) (39b) 

C2 = ((|/3p-^ + 2(Re(a*7))))-^ (39c) 

C3 = (( 2Re(/3*7) - ^Re(a*/3) )) (39d) 

The following cases may then occur: 

• (( |a|2 )) / 0, then co / and dS/clN = 4. 

• ((lap)) = 0, then cq = ci = (since ((a)) = 0). If in addition ((|/3p)) 7^ 1/4, then 
C2 / and dS/diV = 2. 

• (dap)) = and ((1/3^)) = 1/4, then /3 = /3 and cq = ci = C2 = C3 = 0. Hence 
d5/diV = 0. 

The threshold of separability, Eq. (j20p . can be reached by the end of inflation only in the 
first case. 

Let us consider a second Ansatz where the conservation of ( on superhorizon scales is 
violated by a logarithmic term, 



(40) 



where ((a)) =0. The covariances are given by 

Vcc = IQT {1 + x2 + (( lap )) ln2(x)} , (41a) 

V^^ = -H\C'g\' {x' + (( |«P )) ln(x)} , (41b) 

V^^ = hX\'W + {{\»\'))} , (41c) 
and the covariance matrix has the determinant 



4 V 9 / \\ Q J 

The entropy grows with the rate dS/dN ~ 6. Higher rates necessitate a stronger violation 
of the constancy of Q. This behaviour will be observed below in multifield infiation. 



IV. COUPLING TO ISOCURVATURE PERTURBATIONS DURING 

MULTIFIELD INFLATION 

In multifield inflation, adiabatic and isocurvature linear perturbations can be coupled 
on scales larger than the Hubble radius. As a result, even after their decay, isocurvature 
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perturbations have affected the primordial curvature spectra in an irreversible way. In 
this case, tracing over the isocurvature perturbations furnishes a non zero entropy at tree 
level. This source of entropy was considered in We generalize the analysis to a wider 
class of models while using the method of Section IIIII We clarify the conditions leading 
to a linear growth of the entropy with the number of efolds. In particular we consider 
the (non-intuitive) limiting case where the curvature of the background trajectory in field 
space is very weak. 



The model 



This subsection contains review material and may be skipped by the learned reader. 
We consider the following class of two-field models [TTJ 



(i^xyj—g 



R 



levrG 



(43) 



with a non standard kinetic term for the field X; thereby generalizing the action considered 
in The Klein-Gordon equations for the homogeneous background fields are 



26 T 



0. 



(44) 



where the subscripts designate a partial derivative with respect to the field indicated. 
The Einstein equations are 



8ttG 



a 



+ V 



(45) 



where we introduced the field 



■ 2 _ • 2 I 2b -2 

a=ip^ + e^°x ■ 



(46) 



The linear cosmological perturbations are perhaps most transparently written in the in- 
stantaneous basis 5a and 6s of perturbations respectively tangent and orthogonal to the 
background trajectory in the field space ((/?, x) uM 



5a = cos{6)5ip + sm{6)e^dx , 5s = — sm{6)5(p + cos{6)e''5x , 



(47) 



where 



cos ( 



a 



smf = 



e'x 



(48) 



5a is the adiabatic component of the vector of linear perturbations, and 5s is called the 
entropy component. In the present class of models, the anisotropic stress vanishes and the 
line element in the longitudinal gauge simplifies to 



ds' 



(1 + 2$) dt^ + a^{l- 2^>) Sijdx'dx^ 



(49) 



where $ is the gravitational potential. The gauge invariant curvature perturbation (in the 
comoving gauge) is 



H 



$ H 5a , 

a 



(50) 
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where the second expression is obtained using the perturbed momentum constraint equa- 
tion $ + = 4TrGa6a. One also introduces the dimensionless isocurvature perturbation 
during inflation 

S = ^5s, (51) 
a 

which is gauge invariant by construction. 

The linearized equation for (" is remarkably simple, 



where e = —H/H^. The equation for the isocurvature perturbations is a little bit more 
involved 



" dln(a'^e) ■ 
at 

The time dependent function 



(?2 

+ Css 



g \2 $ 



2 



Css = Vss + - b,^^d^ + b%g{t) + b,^f{t) + ^H-^ + ^i-^(^i 
fit) = V^il + sin^e) -AVssine, 

g(t) = -a^ (1 + 3sin2 0) , (54) 
depends only on background quantities. Finally, and S are coupled via the function 

9 ^ -2;^ = |:(^ + &>sm(0) 

~ 2ris^-2b^^sm^i9). (55) 
11 

The last expression is valid in the slow-roll approximation, and rjsa is the slope parameter 
of the potential 

= ^ • (56) 
In brief, C and S are correlated when the background trajectory is curved (in field space). 



B. Qualitative discussion 

Equations (j52p and (|53p are sufficient to understand qualitatively the possible be- 
haviours of the entropy. These equations differ in two qualitative ways, through their 
effective mass and the source. Let us first consider the r.h.s. of the equations. The source 
^ of the adiabatic perturbations is a linear combination of S and its time derivative, but 
does not contain spatial derivatives. It cannot in general be neglected on superhorizon 
scales. The source of S on the contrary is second order in the gradiant of <I> and decays 
exponentially fast. It can therefore be neglected in the long-wavelength approximation. 

Let us now examine the homogeneous equations. While C is massless, and therefore 
reaches a constant value outside the horizon, isocurvature perturbations tend to decay 
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because of the effective mass More precisely, on superhorizon scales (in the slow roll 
approximation), equation ([53|) gives 



T's = Vli-] -I- , (57) 



^HJ \aHJ 
where 

For C55 > 9i/2/4, decreases like 1/a^ which is the damping factor of a massive field. 

The entropy gained by ^ after tracing over S follows therefore one of two behaviours. 
If Css ^ 9H^/A, the isocurvature component remains unexcited (no parametric amplifi- 
cation). The entropy gain is therefore negligeable, frozen at its value at horizon exit, 

We have not attempted to calculate this value. In general, is depends on the integrated 



effect of e, b^^^, and the slope parameters rjij. The analysis of the last reference in [17 | 
shows that in a wide class of models, the slow roll approximation yields good estimates 
of the power spectra at horizon exit. In this class, they depend only on the values of the 
parameters at horizon exit. Unless some fine tuning of these parameters, we expect the 
amplitude of 5* to be 

S. = 0{vl,e^b%). (60) 

In this case, the reduced state of C remains quantum mechanically entangled since 5* <C 
5scp, see Eq. ([20]). 

The other case where C55 < 9H^ /4 is much more interesting. Below we now show that 
the entropy grows linearly with the number of efolds as long as the product gS decays 
slower than 



C. Entropy growth when isocurvature modes are excited 

The analysis simplifies considerably in the long wavelength limit as the damping kernel 
is negligeable. Indeed, in order to get the equation (p2]) for this model, we write the 
solution of (|53|) as 

S{t) = S,{t) - 2 l'jt'Gidt,t') (v,s ^) , (61) 

is the retarded Green function of ([53|) . Sh is solution of the homogeneous equation, 
fixed by choosing the adiabatic vacuum for q/aH :$> 1. Substituting this solution into 
([52]) . one has 

C-^c4c-Al(^/Vc5,„0(,.(^)'i)).«.) (0. 

where ^ is given in (|52p with the substitution S Sh- By identification with (|22p . we see 
that D is suppressed by a factor (^)^. We neglect it in the following with the resulting 
simplification that (j62p reduces to (j52p . Therefore all the expressions below are only valid 
for q/aH < 1. 
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The solution of (j52p has the following form. Before horizon exit, it is given by Eq. ([6]). 
Then it is given by 

m = r{t) + £dt' Gret(t, mt') , (63) 

where C™(^) is given by ^ and G^ct is the retarded Green function of (l3|). To evaluate 
the integral of (|63p . it is sufficient to use the long wavelength approximation ()B6P of the 
retarded Green function. More detailed expressions are given in Appendix IB 21 

To estimate qualitatively the impact of the source ^, we first consider a simplified model 
where ^ is constant during e- folds after horizon exit and where it vanishes afterwards. 
This assumption may appear unrealistic at first, but in the light of the analytical solution 
of the next section, we will see that it contains the essential physics to understand the 
evolution of the entropy. We therefore consider the source term 

^{t) = em-t)eit-Qc, (64) 

where the constant value ^ is given by 

Cc^SH^gSo. (65) 
Then, for < t < t^, C of ([631) tends to 



C(t) = C{t) + gSoln 



a 

at) = C'"(t) + gHSo . (66) 



(Notice that this could have been directly derived from the evolution equation 

+ gS (67) 



dN KaHJ 



assuming S = cte = Sq and neglecting of the first term.) 

Taking into account that at horizon exit, the power spectra are 



2 



„n „n 4:TtG 

(68) 

the covariances for > t ^ t^, are given by 

= Vl [1 + g'^ In^ (x)] , (69a) 

V^^ = + g^H'^Vl . (69c) 
Recalling the relation ([7]) between vr and C, we get 

det(C)c.i+52(^^y [l + 0{g'x^g'ln^x))] . (70) 

We left the subdominant 1/4 to remind that det(C) > 1/4 by Heisenberg uncertainty 
relations. Hence, within the interval 

HU + I In 4^2 1 < m < Ht^ , (71) 
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the entropy grows with a rate 

This expression is rather remarkable since it is independent of the parameters of the model 
and of the initial conditions of (p and x- These appear in the expression of the entropy 
only as logarithmic additive constant, e.g. 21n(g(), and in the boundary of the domain of 
(j7T]) . These features seem to be a generic property of classically unstable systems. 

After the decoupling of curvature and isocurvature perturbations, we verify that the 
entropy is again constant. The expressions are given in Appendix iBl 



D. Canonical kinetic term 

To confirm the physical relevance of the above result, we study in more details the 
dependence of the entropy on g and on the properties of the background trajectory. To 
this end, we turn to the case 

^^ = 0, V = \[mlx^ + m'y) , (73) 



for which analytical solutions exist in the slow roll approximation [19(]. We shall also ex- 
plain the counter-intuitive result of [1], namely that in the limit of a small mass difference, 
the entropy grows linearly till the end of inflation with the rate (j72p while the background 
trajectory is almost straight. 

In the slow roll regime, the background trajectory can be written in the parametric 
form 

X = 2Mpi^sina, ip = 2Mpl^/s cos a , 0<a<^, (74) 

where 



s 



Forward time propagation corresponds to decreasing values of s. The slow roll regime 

corresponds to s » 1 and ends at s ~ 1. We call x the heavy field and introduce the 
relative mass difference 

2 2 

mt — m „ 

X = ^^—^>0. (76) 



The evolution equation 



da A sin2Q ^^^^ 



d In s 4 1 + A sin^ a ' 
is obtained from the Klein-Gordon equations (j44p . Its solution can be written 



.2 



where t = tan(Q;). The index refers to the initial conditions. Prom the solution written 
in this form we immediately see that a is a monotonously growing function of s. Omitting 
the contribution of (j^/2, the Hubble parameter is in turn given by 

H^{s) = \m\s (1 + A sin^ a) . (79) 
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The coupling parameter g has the form 



d6 Xml ^ ^ ^ 



where 9 and a are related by 



1 + 2 



da 1 
d In s tan a 



tan0 = tana din^tani^ , (81) 

l-23t^tana 

Substituting ([77]) into §1^) yields 

tan 6* = (1 + A) tan a, (82) 

which is exact in the slow-roll approximation. In the limit of small mass difference A <C 
1, the velocity vector is almost aligned with the position vector, hence the background 
trajectory is almost straight in this case. 

After horizon exit, the isocurvature perturbations are given by 

where Sq is as previously the amplitude at horizon crossing given in the first approximation 
by (|68p . The gravitational potential is determined by the constraint <^ + = AirGada. 
Together with (jSOp it implies that still reaches a constant and that 

|:^C(l+0(e)) . (84) 



C-gS. (85) 



Inserting this solution into (j67j) . one gets 

ds \aH) 



The first term on the r.h.s. can therefore be neglected. Since gS decays like a power of s, 
which is integrable at s = 0, the amplitude C,{s) indeed asymptotes to a constant. In the 
limit A ^ 1, integrating equation ([55]) yields 

rd(s) 

as) -Co = SoJ^{s,s,) = -2So / deTs{s{9),s,) , (86) 

Je, 

where Co is the amplitude at horizon exit and we used (jSOp and (|83p . This integral is 
manifestly convergent but its expression will not be needed here. Therefore the three 
covariances are given by 

Vcc = V'^ + + J'\s)) , (87a) 

r^i=r^^^ + ro{grsf . (87c) 

V'^^ and V^^ are given by Eqs. (jl2b ) and (c). The part of C driven by S decays polynomially 

in s whereas the homogeneous solution decays as (s/s*)^ exp(2(s — s*)). In consequence, 
in the determinant the term Vq {gTs)^ rapidely overtakes the others, and we have 

1 



det(C) _ ^ 
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where gTs depends polynomially on s. Hence as in the previous section, we find that the 
entropy depends only logarithmically on parameters and initial conditions of the theory 
via In (7* and Insg- In other words, g needs to be exponentially small in order to suppress 
the growth of the factor {aH/qf' . 

More puzzling is the counter-intuitive result [H] implying that the smaller A, the smaller 
the coupling g but the higher the final value of the entropy because the longer it grows. 
We can understand this qualitatively in the following way. Examination of (j54p shows that 
in the slow-roll approximation, the leading contribution to the mass of S is the curvature 
of the potential, 

Vss = ml{l + \ cos2 e) < ml (1 + A) . (89) 

The Hubble rate is such that ZH^ /2 > mis, where the slow-roll condition requires s ^ 1. 
Combining these two inequalities, we get 

9H^ ~ 9H^ - 3 s ■ ^ ^ 

Recalling that isocurvature perturbations decay slowlier than when 4:Css/9H^ < 1, 
we see that the smaller A, the easier it is to satisfy this condition. For A sufficiently small 
it may be verified until the end of the slow-roll regime. On the contrary, for A ^ 1, 5 
ceases to be a source of entropy 0(A) efolds before the end of inflation. In the opposite 
regime of extremely weak coupling, A ^ 1, 5 oc A only enters in a logarithmic shift of the 
time lapse where ([72]) is valid, see ([7T|) . 



E. Discussion 



A high entropy growth rate and the logarithmic dependence on the parameters implied 
hy S = Indet(C) and (j88p seem to be generic for systems whose equations of motion 
are characterized by an instability. For cosmological perturbations, the instability is the 
mechanism of parametric amplification itself. Recall indeed that = x a^, where 
each factor is brought by the relation ([7]) between vr and (. When ( is decoupled, two 
factors of 1/a'^ coming from the Wronskian of C cancel this thereby giving det(C) = 
cte. Instead, in the presence of a source ^, the change in det(C) is given by times 
{T'^5V^^ + V'^-dV(^ - 2V^^^5V (-(-). This product rapidely asymptotes to a^V^dV^^ unless the 

source decays faster than In other words, the conditions of an efficient decoherence 

are provided by the same mechanism producing the highly non-classical state in the first 
place. 

When the growth of entropy follows Eq. (j72p till the end of inflation, the quantum 
coherence of the reduced state is lost since the resulting value of the entropy is 5cnd ~ 
6 In A'end which is higher than the threshold of separability given in Eq. (j20p . We recall that 
this large entropy is however compatible with the classical coherence of the distribution 
(which implies the presence of the acoustic peaks) because it is much smaller than the 
thermal entropy = 8 In A'cnd which characterizes the incoherent distribution (with no peak) 
0- 

In the model studied in Sec. IIVDI the law ()72p is 'generic' in that lower rates < 
dS/dN < 6 are only found when ACss/^H"^ approaches one before the end of inflation. 
Once ACss/'^H'^ > 1, the entropy saturates at 5 ~ 6N^, with 5* = 0(77^^, e^, 6^^) 
and is the number of e- folds from horizon exit till that moment. Hence the threshold 
of separability is reached when > lA'cnd- From the inequality (j90p . this translates into 
an upper bound on the relative mass difference A. To get a rough estimate of this bound. 
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we simplify ACgs/dH'^ by its upper bound ([90|) . The condition 2(1 + A)/3sg = 1 with 
= Nend — ^s, then gives the upper bound 

^ e-3^« < A < ^ - 1 . (91) 
2 - - 2 ^ ^ 

The lower bound on A comes from (f7T]l and ^ A/s* = X/Ncnd- It is easily satisfied 
when the number of e-folds N^^ is not fine tuned (e.g. by the choice of initial conditions) 
to be small. 

It is simple to extend our conclusion to inflationary models with M fields. Then, ^ is 



the sum of the contributions of the M — I isocurvature modes |2ll |. Hence, the entropy 
grows steadily as long as at least one isocurvature mode is excited. 

It is harder to reach conclusions in the case of multifield inflation with a non-canonical 
term. In this case, the derivatives of b{Lp) can contribute significantly to C55 so as to 
increase or lower the effective mass. Moreover, the coupling can be strong enough to 
lead to a significant change of C on superhorizon scales and to invalidate a perturbative 
treatment. A numerical treatment is probably required to settle the question. 



V. COUPLING TO MATTER FIELDS 

In a non stationary background, there are a priori two distinct sources of entropy from 
loop corrections. The first is dissipation, as in non vacuum states in Minkowski space. 
The second is a possible non trivial time dependence of radiative corrections that would 
frustrate the cancellation between the variances which at tree level lead to det(C) = 1/4. 
Weinberg identified the conditions such that the corrections depend only on the values 
of background quantities at horizon crossing (instead of the whole history after horizon 
crossing). These conditions are satisfied by a wide class of semi-realistic theories 
We verify that in these theories, the entropy vanishes at one loop approximation. 



The model 



We consider theories of one inflaton field (p{t,x). Matter is modeled by copies of 
free scalar field a. More precisely, we consider the vector field a{t,x) = {ai, ...,aj^) in 
the fundamental representation of 0{M). We assume no symmetry breaking otherwise 
isocurvature perturbation are excited. The gravitational sector was described in Sec. HH 
The gravitons have a similar action but as explained below we will not need to take them 
into account. The quadratic part of the action describing the free evolution of the matter 
fields is 



S[crn] 



1 ^ r 

n=l 



dtcfx 



■2 



1 



-o (V(T„) 



(92) 



Minimal (conformal) coupling to gravity corresponds to = (resp. ^ = 1/6). 

C and the matter fields are canonically quantized. We work in the interacting picture. 
The momenta conjugate to (T„ are 



vri(t, 



,3 ■ 



0-j(t, 



(93) 



and the mode equations are 



d\n{a^ 
dT 



-<yq + 



(94) 
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The free vacuum is the Bunch-Davis vacuum defined by 

q 

(idr — q) (ao-q) — > for > oo (95) 

aH 

In a quasi-de Sitter approximation, the mode functions are 

= aO-''^(l + zgr)e-^''^ (96) 
ag{t) = aO-°fe-^\ (97) 

for minimally and conformaly coupled scalars where the normalization constants are 



^0 mini 2 _ q |_0conf|2 _ /Qn\ 

I Wq I -Yq- (98) 



We work at the leading order of the large- limit. At this order, inspection of the dia- 
grammatic expansion reveals the following elements. First, gravitational self-interactions 
are irrelevant. It means that the gravitons decouple from the scalar perturbations, and 
that scalar self-interactions do not contribute. Second, since each matter loop is enhanced 
by a factor M, matter fields propagate freely. Therefore the only loop corrections to con- 
sider are matter loops in the two-point function of C- Moreover, only the trilinear vertex 
(aa contributes to the logarithmic part of the one loop correction (the local-regular part of 
the loop correction does not contribute to the entropy [ij], so that the one-loop correction 
comming from the vertex CC^'"'^ need not be considered here). Finally, an explicit calcula- 
tion 0] shows that the relevant part of the trilinear vertex responsible for the logarithm 
is 

S^„„ = - jdtHi^t = - J dtd?x (r°° + o^SijT^) (^-eHa'^V-'^c) , (99) 

where T^j/ is the energy momentum tensor of matter. For minimally and conformaly 
coupled scalars, the linear combination in (j99p is respectively 



(r00 + a25,,r^)"^'" = 2&\ (100) 
(TOO + a%r^r''' = cr' + l (^(Va)2 - 2aa - H^a'^ . (101) 



B. Outline of the calculation 

The expectation value of a local (possibly composite) operator Q{t,x) is given by 
Q(t,x) = Q/(t,x) +i I dt2[Hj{t2), Qi{t,x)] 

J —oo 

- [ dt2 r dti [Hj{ti),[Hj{t2), Q/(t,x)]] + ... (102) 

J —oo J —oo 

The subscript / refers to the interaction picture and will be omitted in the following. In 
this expression, Hj is interacting Hamiltonian, e.g. H[^t in (|99|) . The dots stand for higher 
order corrections and counterterms (whose explicit form will not be needed here). The 
Fourier transform of the one loop correction to covariances SV^^i can then be written 

ld'xe^^-5m^)^'it,0)) = -4Aa|^|^ (2^)35(3) (q+p + p') 

X / dt2V2 dtiViRe(zfM^p,) + ... (103) 

J —oo J —oo ^ ^ 
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where ^ and ^' can be either C of C, and V{t) = -eHa^. The same function Ai appears in 
the calculation of the three covariances. It depends only on the matter fields an and the 
loop variables p, p', ti and t2- If the cr's are conformally coupled, one finds 



M 



conf 
PP' 



M 1 



p2 _|_ ^' 2 



App' 



36 a\a\ pp' 
For minimally coupled scalar fields one gets 



-i{p+p')(Tl-T2) 



{pp') 

4 4 



-i{p+p')(n-T2) 



(104) 



(105) 



Notice that the time dependence of these two expressions is the same. This remark is 
essential to understand why the entropy from minimally coupled fields is not much larger 
as might have been naively expected (since they are not in the conformal vacuum but 
are parametrically amplified). The function Zq^ depends only of the modes C,q and/or 
Using the solution ([6]) of the mode equation, we get 



a^{t)H HiH2alal 

ircc = __S. I'-gi 



g-ii3(ri+T2-2r) _ g— «(j(ti-T2) 



(106a) 
(106b) 
(106c) 



In each line, the first term in the brackets comes from the term (C1C2Q) of the perturbative 
expansion (jl02p . and the second from (CiQC2)- It is tempting to dismiss the terms 0{qT) 
since for the power spectrum of super Hubble scales, they represent subleading terms. 
These must however be kept to calculate the entropy. Indeed, we remarked already below 
Eq. (fT5]) that at tree level, three powers of 1/gr cancel in the expression of det(C) so as 
to ensure that the entropy vanishes, see Sec. HIl 

Then, the key observation is that all the factors of a coming from V{ti^2)^ ^ a-nd M 
cancel in the integrand of (jl02p . Thus, the dummy variables ri and T2 appear only in the 
phase of exponentials. Explicitely, they are the phase factors in the brackets of Eqs. p06|) 
and from M of Eqs. (jl04p and (jlOSp . The singular logarithm of the one loop correction 
can be calculated by exchanging the order of integration. The integration of the first phase 
in Eqs. (jl06p gives 



ii{q,p,p') = e 



,i2qT 



/T2 
-00 



dri e 



-i{p+p' +q)Ti 



2q{q + p + p') 



.(107) 



The remarkable property of this term is its r-independence. This is a consequence of the 
stationarity of de Sitter space, but it can be shown that the integral is also finite in the 
limit t +00 in power law inflation for e < 1/3. For the second term, we introduce 
Q = p + p' + q, and make the change of variables {t2,ti) i— > {t2,St = ti — T2). We get 



'^2{q,p,p') 



T Q + ie 



At 



(108) 



where At = t + T. 

The integrals I\ and I2 times the kernel M must then be integrated. Since > 0, 

the Dirac 6 in (jlOSp gives a vanishing contribution for any finite while the principal 
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value gives a purely imaginary contribution. Since this second term is only multiplied by 
real functions in Eqs. (jlOOp . it therefore does not contribute to (|103p . Hence, only the 
interfering term {C1C2Q) contributes. Because of the the curvature perturbation couples to 
matter via its time derivative C, its conjugate momentum is not d?]) but vr = a^e(^/(47rG) + 
F{da, &) where F is a functional quadratic in a. These additional terms do not contribute 
at one loop to the logarithmic part of the variance "P^^ and "Ptttt so that we ignore them. 
Therefore the 1-loop modifications of the variances are 

(Cc(t)C-q(t))i-loop =M {Se'X/) (1 - ^) J (^) , (109a) 
liiUt), 7r„q(t)})i_ioop = -M {seXf) (--^ ^) J , (109b) 

(^,(t),vr_,(t))i„,oop=AA(86^|C°r) (-J^^)'^(f) ' (109c) 

where J contains the (ultra-violet divergent) integral over loop momenta. For instance, 
for a minimally coupled field its un-subtracted expression is 

- - /7?T3 7f4(2-)^^^^Uq + P + P') — (110) 

The new scale // in Eq. (jl09p is an a priori arbitrary scale associated with regularization 
procedure, e.g. dimensional regularization. The occurrence of g > in (|107p regulates the 
momentum integral in the infrared. This is the advantage of having inverted the order 
of the integrals. Note that the expressions (jl09p hold for both minimally and conformaly 
coupled fields. We call A the coefficient of the logarithm of J . Its value depends of course 
on whether a is minimally or conformally coupled. It is given in d, 23] but we shall not 
need its explicit expression here. Taking into account the normalization ([8]), the relative 
change of the power spectrum is 



5V^ 



ANeGH'^ln , (111) 



up to terms (q/aH)'^. 

The covariance matrix does not grow with a certain power of a as it could have done 
a priori. It is given by 

det(C) = \ + {egGH^MA ln{q/fi)) ^ . (112) 

In conclusion at one loop approximation, 

S = 0{eg X 10-^° X AA)2 , (113) 

even though the power spectrum is modified to linear order in eg x GH^ x M. Comparing 
with the entropy of a separable state in the radiation dominated era Sscp = lOOln(lO), the 
reduced density matrices remain entangled (not separable) at the end of inflation, unless 
the number of fields is larger than l/e^ x lO^'' but in this case the whole perturbative 
treatment is no longer valid. 

The reduced density matrix of modes of opposite wave vectors remains very pure, and 
can thus be interpreted as a "dressed squeezed vacuum state". This is an example of 
the case discussed in Sec. IIIICI We stress that this conclusion could not be expected 
a priori on the basis of an analogy with field theories in the Minkowski vacuum. In 
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inflation, matter fluctuations of non conformal fields are parametrically amplified just 
as the curvature perturbations, so that they could have been responsible for a strong 
decoherence. As we saw, this is not the case and mainly follows from the fact that matter 
couples to C via its energy-momentum tensor, see ()99p . 



VI. CONCLUSION AND OUTLOOK 

We found that in multi-field inflation scenarios, the entropy grows at a high rate, 
typically dS/dN = 6, as long as some isocurvature mode decay polyniomaly in the number 
of efolds after horizon exit. By contrast, in single field inflation we found (at one-loop) no 
evidence for decoherence, so that the state remains essentially pure. Intermediate rates 
between these two cases seem hard to find because of the efficiency of the mechanism of 
parametric amplification (in amplifying or canceling the action of a noise). 

We concentrated on the entropy of curvature perturbations, but our method may as 
well be applied to tensor modes. We conclude on the perspectives to apply this method 
to the regimes of (p)reheating and horizon reentry. One should distinguish between sub- 
and superhorizon scales. The evolution of the former is described by a regime of broad 
parametric resonance. This is a process characterized by its efficiency. Decoherence should 
therefore be achieved very rapidely. The evolution of super-Hubble modes depends on 
whether isocurvature perturbations are excited or not. The discussion of this case can be 



incorporated in the model of two field inflation of Sec IIVI (see 2J] as well as references 
therein). 

After reheating but before decoupling, it is expected that the coupling of curvature 
modes entering the horizon with the radiation density perturbations in quasi-thermal 
equilibrium will erase the possible remaining quantum features after typically one oscil- 
lation (this should not be mistaken with the thermalization of the perturbations which 
occurs much later). We have little to say about this regime of horizon reentry. A field 
theoretic proof remains a formidable task because long wavelength radiation perturbations 
are not thermal fluctuations since we observe acoustic oscillations. Since this phenomenon 
is reliably described by a hydrodynamic model, it means that scales decouple. As a result, 
we expect that the curvature perturbations entering the horizon dominantly couple to the 
same scales of the plasma. A proper investigation of the regime of horizon reentry there- 
fore requires first to write an effective field theoretic model of these out-of-equilibrium 
long wavelength modes. This is arguably an academic question, since we saw in [H that 
the loss of quantum coherence occurs on time scales much shorter than the characteristic 
time of thermalization. 
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APPENDIX A: DERIVATION OF EQ. ([22]) 
1. Equation for the anti-commutator 

When the Hamiltonian depends exphcitly on time, there is no stable ground state. In 
this case, Green functions, e.g. 

G{x, y) = (0in|rv9(x)v9(y)|0in) , (Al) 

are expectation values in the 'in'-vacuum |Oin) which cannot be expressed in terms of 
Feynman graphs with internal lines corresponding to time-ordered propagators. Rather, 
their generating functional is given by the transition amplitude of two 'in'-vacua in the 



presence of external sources J+ and J_ (see for instance j25l |) 



where T is the reversed-time ordered product and J+ and J_ are the two classical sources 
associated with the two branches of evolution, forward and backward in time respectively. 
Note that the operation of taking the trace couples the forward and backward time evolu- 
tions. As a result, 1^[J_|_, J_] generates four types of connected two-point functions: the 
time-ordered propagator 

5W 

G++{x,y) = i{Tip{x)ip{y)) = \j+=j_=o, (A3) 

the reverse times ordered propagator 

5W 

G..{x,y) ^ ^{T^{xMy)) = Jj^^^^^^^W=J-=o , (A4) 
and the two on-shell two-point functions 

6W 

G-+{x,y) =i{ip{x)ip{y)) = —————-\j j^=o, (A5) 

6J^{x)5J+{y) 

6W 

G+^{x,y) = i{ip{y)ip{x)) = A j+=j_=o- (A6) 

6J+{x)5J^{y) 

The latters are not time ordered because they are build from operators coming from 
different branches of time evolution. Eqs. (|A3IIA6p are easily derived using the second 
equality in ()A2p . 

The path integral representation of (IA2p is 



where S[lp'^, (/9^]((/?)^) is the classical action evaluated for the paths (j)~^{z^,'z) with fixed 
end-points (oo,x) = 99+ and (/)+(S,y) = 93+. 

In the above equations, 93 denotes the collection of fields. In cosmological settings, we 
decompose these fields into the adiabatic perturbations C, and the rest we call a. Then, the 
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so-called influence-functional (IF) encodes the effective dynamics of C when these extra 
fields have been integrated: 

/+00 r+co 
-oo J — oo 



We can evaluate perturbatively the IF by expanding the terms e*('^'"*[*+l~'^'"'[*-l) to a 
given order in ^, then by taking the expectation value for the field a, and finally by 
re-exponentiating the result. Hence we get 

5,^[C-,C+] = i{{s,nt[C+])-{suC-])) 

— - { {Sint[C,+\Sint[C,+\) con + {Sint[C-]^int[Q-]) con 

— {Sint[C+\Sint[Q-]) con — {Sint[C-]^int[C+\) con} + ••• (A9) 

where ( )con means the (connected part of) expectation value in pj^, the initial state of 
the field a. 

The Gaussian approximation consists in keeping only the part quadratic in C,± of this 
functional. Then the generating functional is given by 

giW^gau..[J+,J-;Pi„] = jvC^VC exp{-^^Z MZ + ^JZ} 

cx exp (^-^'JM-^J^ , (AlO) 
where *Z = ((^-|_ ^_) and ''J^ = i (J+ , — J_). The quadratic form Ai is by definition 



'ZMZ = 5[C+] - S[C-] + Sfr^iU, C-] • (All) 



1 

2 

The action S is given by eq. ([2]) whereas the Gaussian approximation of the IF is 



Sfr" = ^ jdtd\{a\t)e) Jdt'd^y{a'{t')e) Ca{x)^ab{x, y)Cb{y) , (A12) 

where the indices a, 6 = it. This equation defines the self-energy matrix T,ab- We have 
introduced twice the factor a^e, which it is already present in the kinetic action ([2]), because 
it simplifies the forthcoming equations. 

Taking the functional derivatives with respect to the sources, we arrive at a system of 
two linear coupled equations. The system decouples when using the (odd) commutator 

(also called the spectral function) p = i{G |_ — G-| ) and the (even) anticommutator 

G = (G+++G )/2 = (G-i +G Indeed, using the definition of the free propagators 

= 'Dx = Gq one gets 

Vxp{x,y)+ d^zD{x,z)p{z,y) =0, (A13) 

JyO 

P,G(x, y) + r d^z D{x, z)Giz, y) = f d^z N{x, z)p{z, y) . (A14) 

Straightforward algebra gives the odd and the even kernels 

D{x,x') = ia^{t')[E^+{x,x') -T,+^{x,x')] , 

N{x,x') = ^^[S.+ (x,x') + S+_(x,x')] , (A15) 
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in terms of T,ab, the self-energy matrix of 

We conclude this part by a comment. In the Gaussian approximation it is always 
possible to write the effective action with only the field (here C) but no field derivatives. 
Should one find, after computation of the effective action, field derivatives (e.g. because 
in the total action the curvature perturbation couples to other fields via Q, one would 
simply do an integration by parts in order to bring the effective action into the form ()A12p . 
This operation is a linear canonical transformation which changes neither the equations 
()A13l IA14|) nor the value of the determinant of the covariance matrix Moreover, 
with the action written in the form (|A12p . the canonical momentum is still given by 
vr = (a^e/47rG) C which justifies the identities (fTTj) . 



2. The effective quantum source 



Since the unknown of Eqs. (|A13|) and (IA14[) are real functions, they are apt for a 
numerical analysis. Alternately, for an analytical treatment it is simpler to go one step 
backwards and consider the operator <^q coupled to a quantum mechanical source whose 
statistical properties are such that the anti-commutator of (q obeys by construction Eq. 
()A14p . This effective source is not just an artificial trick, because it coincides with the 
true fluctuating source operator when non-Gaussianities are neglected, as it is the case for 
the two models considered in the body of the paper. 

Let us consider the following Heisenberg equation 



A$q(t)+ / dsD{t,s,q)<!>, 

J — oo 



(A16) 



where and .^q are two operators. We also introduce the anti-commutator 

gg{t,t')6'iq-c^) = ^Tt [p^p^ {c|.q(t), a>t_^,(t')}] ■ (A17) 

For an appropriately chosen anti-commutator of ^q, we now show that given some initial 
conditions, the function Qq{t,t') solves for the (Fourier transform in space) of ()A14|) with 
the same initial conditions. Hence it is equal to the anti-commutator G{t,t',q). 

To determine the statistical properties of the noise ^q(t) that reproduce the action of 
the r.h.s. of jATH), we let the integro-differential operator on the l.h.s. of (|A14p act on 
Qq, and calculate the result using ()A16p . Since the action of taking the trace commutes 
with the partial derivation and integration, we have 

ft 



<5^(q-q': 



VtGqit,t')+ [ dsD{t,s,q)Gq{t,t') ={Uq{t),<^l^,{t')\). (A18) 

J —oo J ^ ^ 

To calculate the anticummutator on the r.h.s., we need the solution of (IA16p . Remember- 
ing that the exact retarded Green function of $q is the spectral function multiplied by a 
theta function, i.e. 6{t — t')pq{t,t'), the general solution of (1A16P is 

^q(i) = ^q(*) + [ ds pq{t, s) (Al9) 



where ^'q(t) is the solution of the homogeneous equation. Since it is independent of ^(t), 
upon substitution of this solution into the r.h.s. of (jAlSp one finds. 



53(q-qO 



Vtgq{t,t')+ I dsD{t,s,q)gq{t,t') 

J —CO 



ds{{^q{t), C-^,{s)})pq{s,t') 



(A20) 



25 



Identification with (|A14p finally yields 

^({^qW, = N{t,s,q)6H^-^) (A21) 

where N{t,s,q) is the (Fourier transform of the) kernel appearing in Eq. ()A14p . In 
conclusion, the dynamics of ^ is equivalent to the "open dynamics" (jA16p of <I> with 
the same spectral function p and subject to the Gaussian source characterized by the 
spectrum ()A2ip . 

This result is general because we are only considering the two-point functions. Because 
of the linearity of the equations and the Gaussianity of the noise, this approach can be 
also phrased in terms of a stochastic (commutating) c-number source ^q, see the stochastic 
approach of quantum gravity in [2^. This equivalence is used in Sec. IIIIDi 

Finally, from the definitions (jlip and the identification G = Q, the covariances can 
always be written as a the sum of a free (and decaying) part and a driven part: 

($°(i)$°(t)) + jdtidt2G,,t{tM)Grct{tM)N{hM), (A22a) 

\{[^l{t), ^l{t)]) + Jdhdt2Grct{t,h)dtGret{t,t2)N{h,t2), (A22b) 

(60(t)^.0(t)) + [dtidt2dtG,,t{t,h)dtG,,t{t,t2)N{h,t2) . (A22c) 



APPENDIX B: INTERMEDIATE RESULTS FOR SEC. |IV] 
1. The free retarded propagator in the long wavelength limit 

The free retarded propagator of ([3]) is 



Grct{t,t',q) = 2e{t-t' 



W{t') 



(Bl) 



where (g and Q are the homogeneous growing and decaying solutions of ([52|) . and W{t) 
is their Wronskian. For our calculation, we only need to retain the leading terms of their 
expansion in powers of (q/aH)'^. In the limit g — > 0, the Wronskian solves for the equation 



dt 



whose solution can be conveniently written as 



W{t) = W{U 



are. 



where is the time of horizon crossing q = a^H^. The solution of ([3]) is 



Co(t) = c; 



1+0 



dt' 



a^{t')e{t') 



+ 



(B2) 



(B3) 



(B4) 



The coefficients and Aq of the growing and decaying solutions are related by the Wron- 
skian condition 



i ( CqCq - C.C 



2Im(C°Mg) _ AttG 



a^{t)e{t) a\t)e{t) 



(B5) 
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Combining these results, we get 



Gretit,t',q) = e{t-t')a\t')e{t') 



+ 



9{t - t') 
3H 



a{t) 



0{t-t') 
3H 



(B6) 



To get the second Hne, we substituted the solution (IB4h of the mode equation in the slow- 
roll approximation. We also give the first time derivative of G^ct in that approximation, 



dtGrct{t,t',q) 0{t-t')(^^ 



i\\3 



(B7) 



2. The covariances of Sec. IIV Bl 

We use the solution ([6]) of the mode equation. The retarded Green function is 

Im{Q{t)Qiti)} 



Gret{t,h) = -2e{t-ti) 



^^^j^ /m { (1 - ix){l + ixi)e^(^-^i)} 



We use again the variable 



The noise kernel is 

N{ti,t2) =9H^g^Vo for t^>ti,t2>U 
The covariances calculated from (|A22p are given by 

rc( = Vo{l + x^ + g^fix;x^)} 

V^^ = -Hx^Vo { 1 - 9^ fix; x^)h{x; xg)} 

V^^ = H^x^Vo{l + 9^h\x;x^)} 

and the determinant of the covariance matrix is 



det(C) = {x^ [1 + 9^h\x; x^)] + 9^ {f - hf] 



(B8) 



(B9) 



(BIO) 



(Blla) 
(Bllb) 

(Bllc) 



(B12) 



The functions / and h come from the integral of Grct(i)ii) and dtGrctit,ti) respectively. 
Their expressions are 

f{x;x^) = Im {{1 - ix)e'''^J{x;x^)} = h{x;x^) - X Re{e'''^J{x;x^)} (B13) 
h{x;x^) = /m{e"'J^(x;xg)} (B14) 



and the function J{x\x^) is 



^dxi 



Ji A 



(B15) 
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a. Entropy growth for t < t^^ 

In that case 



J{x; xA = -iC{x) - e""' ( A + 4 + - ) + e"* (2 + i) (B16) 

\x-^ X"^ X J 

where the first term is a logarithm 

C{x) = T— e-*^i = Ei{i) - Ei{ix) = ln{x) + 0{x) (B17) 
Ji xi 

Combining these expressions, we get 

/ = -ln(x)+Ci + 0(x) 

h = -\- ln{x) + C2 + 0(x) (B18) 



from which we obtain the expressions (|69ti69|) as well as 



2 



det(C) = ^{l + 0(x2)} (B19) 

b. Constant entropy after 

We see on (]B12p that the entropy is constant provided 

x^h^ + (/ - hf = Xx^ (B20) 

where A is a constant. For t > t^, the integral J is equal to its value at 

J{x; x^) = cte = J^ = U^le^^ (B21) 

Substituting into the expressions of / and h, we find 

f-h = -xi2e {e^Jg} = -x| J'gl cos(x + 

h = xlm{e"j^} = x\J^\ sm{x + If) (B22) 

The condition (|B20p is realized and we have 

det(C) = i{l + 5Vd'} (B23) 

with \J^\=x^^. 

One also checks that the covariances (fBTTl) verify the identities ^ and ([25]) toghether 
with the free equation of motion. 
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